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I.     INTRODUCTION 

Compared  with  mechanisms  for  the  growth  of  water  waves,  those 
for  wave  decay  have  received  little  attention.   In  wave  forecasting, 
for  example,  there  is  far  more  agreement  on  empirical  wave  growth 
formulas  than  on  decay  formulas.   In  the  case  of  decay  due  to  ambi- 
ent ("external")  underwater  turbulence,  the  relative  neglect  is 
probably  due  to  the  lack  of  both  ocean  turbulence  data  and  a  simple 
theory.   That  of  Phillips  (1959)  is  based  on  the  turbulence  spec- 
trum near  a  free  surface,  which  has  not  been  studied  in  detail. 
Earlier,  more  speculative  work  (Groen,  1954;  Bowden,  1950)  is  also 
closely  tied  to  unknown  characteristics  of  the  near-surface  turbulence 

However,  wave  decay  is  quite  important,  and  the  results  of  the 
simple  but  novel  laboratory  experiment  described  below  provide  some 
insight  into  this  problem.   These  results  are  admittedly  fragmentary: 
the  experiment  was  preliminary  in  nature,  and  the  work  had  to  be 
abandoned  unexpectedly.   Nevertheless,  it  seems  worthwhile  to  pre- 
sent our  data  and  tentative  conclusions  to  guide  others  in  similar 
work . 

II.   APPARATUS 

A  paddle  (or  fan)  generated  surface  waves  which  propagated 
down  a  long  tank,  over  a  zone  of  underwater  turbulence,  to  a  non- 
reflecting  beach  (Fig.  1).   The  turbulence  was  produced  by  a  ver- 
tically oscillating  underwater  grid,  which  spanned  the  width  of  the 
tank,  and  consisted  of  l.-cm  wide  metal  strips  woven  about  1 . -cm 
apart  so  that  the  energy  containing  eddies  had  a  scale  of  1.  cm 
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(see  Paquin,  1968,  for  details).   The  oscillation  frequency  was 
approximately  1.  Hz;  the  amplitude  was  1.  cm.   The  eddies  had 
little  direct  effect  on  the  water  surface.   The  waves  produced  by 
turbulence  alone  were  less  than  a  millimeter  in  height,  whereas 
the  paddle-  or  wind-generated  wave  heights  were  on  the  order  of 
centimeters  upstream  of  (i.e.,  before  entering)  the  turbulent 
region.   No  wave  breaking  was  observed. 

Surface  waves  were  shielded  from  mean  currents  caused  by  the 
grid  drive  mechanism  by  a  false  bottom  (Fig.  1) .   Any  such  current 
remaining  was  completely  masked  by  the  velocity  fluctuations ,  so 
that  the  effect  of  wave-current  interactions  was  probably  quite 
small.   Wave  frequencies  ranged  from  1  to  10  Hz.   Frequencies  above 
2  Hz  corresponded  to  deep-water  waves  (Stoker,  1957),  so  that  the 
influence  of  both  bottom  and  grid  was  negligible.   For  lower  fre- 
quencies the  false  bottom  was  removed  (introducing  weak  currents 
not  otherwise  present)  but  the  waves  did  feel  the  grid  itself. 
Faraday  waves  (McGoldrick,  1968)  were  suppressed  by  thin  metal 
strips  near  the  paddle. 

Resistance  wave  gauges  were  mounted  upstream  and  downstream 
of  the  turbulent  zone,  and  were  calibrated  both  statically  and 
dynamically  (using  progressive  waves  and  stroboscopic  light) .   The 
response  was  flat  from  DC  to  4  Hz  and  10%  down  at  26  Hz.   A  con- 
stant-temperature anemometer  (DISA  Model  55A01)  was  used  with  a 
quartz-protected  hot-film  probe  (DISA  55A81)  to  measure  the  turbu- 
lent motions.   All  data  were  FM  recorded  with  a  Precision  Instruments 
6200  Tape  Recorder  and  analyzed,  after  playback  at  10  times  or  100 
times  the  recording  speed,  with  a  General  Radio  Type  1900  Wave  Analyzer 


III.   MEASUREMENTS 

All  measurements  were  made  during  steady  conditions  with  clean 
tap  water.   The  decay  due  to  viscosity  and  possible  surface  films 
was  determined  first,  then  (immediately  afterward)  that  due  to  these 
causes  plus  turbulence. 

A  typical  paddle-generated  wave  spectrum  (with  or  without  turbu- 
lence) is  shown  in  Fig.  2.   In  order  to  investigate  the  widths  of 
the  harmonics,  the  recorder  playback  speed  and  the  analyzer  bandwidth 
were  varied  to  produce  an  effective  bandwidth  that  ranged  from  1.0  Hz 
to  0.1  Hz.   Since  the  observed  energy  level  did  not  depend  on  the 
bandwidth,  the  wave  energy  was  assumed  to  be  located  only  at-  the 
paddle  frequency  and  its  harmonics,  and  results  were  interpreted  as 
line  spectra.   Figure  3  shows  typical  line  spectra  in  the  presence 
of  turbulence. 

The  turbulence  velocity  spectrum  obtained  with  the  anemometer 
placed  in  the  center  of  the  turbulent  region  at  a  depth  of  2.  cm,  is 
shown  in  Fig.  4.   Similar  spectra  were  found  at  depths  of  1.  and  3. 
cm  at  various  positions  in  the  region.   The  lack  of  a  large  super- 
imposed unidirectional  flow  introduced  serious  problems  associated 
with  variations  of  probe  sensitivity  with  flow  direction.   Calibra- 
tion showed  the  minimum  sensitivity  to  be  about  25%  of  the  maximum. 
(The  hot-film  probe  was  always  set  to  be  most  sensitive  to  flow  in 
the  wave  propagation  direction.)   Thus,  the  turbulence  spectrum 
should  be  interpreted  in  a  qualitative  sense,  showing  only  that  the 
spectral  peak  is  below  the  wave  frequencies,  and  giving  some 
indication  of  the  turbulence  intensity. 


A  corresponding  velocity  spectrum  due  to  paddle  waves 
alone  is  also  shown  in  Fig.  4.   Despite  the  sensitivity  prob- 
lem, the  peaks  associated  with  the  paddle  waves  dominate,  and 
the  energy  associated  with  the  fundamental  wave  is  within  a 
factor  of  two  of  that  predicted  by  linear  theory. 

Within  the  interacting  region  the  turbulence  energy  is 
comparable  to  the  wave  energy.   This  is  convenient  for  the 
experiment  since  large  (hence  easily  measured)  wave  decay 
results.   On  the  other  hand,  it  probably  does  not  reflect  con- 
ditions in  a  natural  environment,  where  wave  motions  are 
expected  to  dominate  (Shonting,  1968). 

IV.   DISCUSSION 

A.   Paddle-generated  Waves 
According  to  laminar  viscous  theory  for  surface  waves 
(Van  Dorn,  1966)  ,  wave  decay  with  distance  is  given  by  the 
linear  expression 


da/dx  =  -B„a 


which  gives 


a(x)  =  a(x  )exp{-B.(x  -  x  )  }, 
o  X,         o 


(1) 


Here,  a  is  wave  amplitude,  x  is  propagation  distance,  and 
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the  damping  coefficient,  varies  with  frequency  f,  viscosity,  and 
the  character  of  any  surface  film  present.   This  is  a  first-order 
result  in  ka  (k  is  wave  number) ,  valid  for  ka<<l.   The  measured 
ka  values  ranged  from  0.04  to  about  1;  most  were  near  0.5. 


Nevertheless,  we  use  (1)  to  determine  8  .  In  no  case  was  the 
observed  laminar  decay  greater  than  could  be  explained  by  the 
presence  of  a  surface  film. 

Following  Phillips  (1959),  two  turbulent  mechanisms  can 
lead  to  wave  decay.   Scattering  (due  to  turbulent  convective 
distortion  of  wave  fronts)  is  a  first-order  process  in  ka , 
results  in  the  generation  of  a  broad  spectrum  of  scattered 
waves,  and  leads  to  an  exponential  decay  law  analogous  to  (1) 
with  3 „  replaced  by  6. •   Phillips'  results  can  be  rewritten  to 
show  that  two-dimensional  waves  passing  through  a  horizontally 
homogeneous  turbulent  region  experience  scattering  losses 
giving  &°<f 

Dissipation  (wher-jby  mean  wave  motions  strain  turbulent 
vortex  filaments,  thus  increasing  turbulent  energy  at  the 
expense  of  wave  energy)  is  a  second-order  process  in  ka ,  which 
is  effective  when  the  scale  of  the  turbulence  is  much  less  than 
that  of  the  waves.   The  decay  law  associated  with  this  mechanism 
is  not  known,  but  is  probably  nonlinear  in  wave  amplitude  (in 
the  differential  form). 

It  is  likely  that  dissipation  dominated  in  this  experiment 
for  the  following  reasons:   (1)  the  mean  square  surface  height 
downstream  of  the  turbulent  zone  was  typically  less  than  10%  of 
that  upstream.   (2)   Any  decay  due  to  scattering  should  be  accom- 
panied by  a  broadening  of  the  downstream  wave  spectrum;  we  could 
not  detect  any  broadening  with  our  filters.   (This  supported  the 
observation  that  the  downstream  waves  continued  to  appear  long- 
crested.)   (3)  Finally,  scattering  should  cause  random  fluctuations 


in  the  contribution  of  any  fixed  frequency  band  to  the  rms  sur- 
face height,  as  read  on  the  wave  analyzer  voltmeter.    (Note  that 
the  tank  walls  were  smooth,  allowing  good  reflection,  so  that 
wall  effects  on  scattering  should  be  slight.)   Fluctuations  of 
about  10%  of  the  mean  value  were  observed  downstream;  they  were 
undetectable  upstream.   Therefore,  scattering  accounted  for  a 
very  small  fraction  of  the  total  energy  lost  by  the  incident  wave 

B.   Wind-generated  Waves 
Wind  waves  generated  in  the  laboratory  obey  an  exponential 
growth  law  (Hidy  and  Plate,  1966) 

da/dx  =  3  a. 
w 

Since  both  the  growth  due  to  wind  and  the  laminar  decay  follow 

an  exponential  law,  they  can  be  combined  to  give  an  effective 

3   in  (1)   (which  may  no  longer  be  positive) . 

The  energy  is  now  broad  band.   Measurements  and  calculations 

were  done  using  energy  in  1.-  Hz  bandwidths .   The  results  are 

plotted  against  the  center  frequency  of  the  filter  band. 

IV.    INTERPRETATION  OF  RESULTS 


The  decay  measurements  are  interpreted  in  two  ways. 

A.   Linear  Decay  Hypothesis 
Figure  5  shows  the  effect  of  assuming  a  linear  turbulent 
decay  process  (such  as  would  occur  with  scattering) ,  which  is 
independent  of  laminar  decay: 


or 


da/dx  =  -( B   +  6  )a 


a(x)  =  a(x  )  exp  {-(60  +  8.)U  -  x  )  }. 
o  it  o 


(2) 


The  observed  decay  with  no  turbulence  gives  3   (eqn.  1) ;  the  cor- 
responding decay  with  turbulence  then  gives  3  . 

The  fundamentals  travel  as  free  waves.   The  second  harmonics 
are  not  free  waves,  but  bound  waves  associated  with  a  periodic 
nonsinusoidal  profile  (Lamb,  Art.  250).   We  are  interested  in  the 
decay  of  free  waves.   A  reasonable  correction  is  to  plot  the 
second-  harmonic  turbulent  damping  coefficients  against  the  fre- 
quency of  free  waves  of  the  same  wave  number.   In  the  same  spirit 
this  correction  is  based  on  linear  theory.    (Aside  from  some 
third-order  work  by  Wilton,  1915,  there  are  no  nonlinear  disper- 
sion results  for  capillary-gravity  waves.)   A  free  wave  obeys  the 
dispersion  law 


4iT2f2    =    gk    +    Tk3/p 


(3) 


where  g  is  gravity,  p  density,  and  T  surface  tension.   Since  the 
second-harmonic  wave  numbers  and  frequencies  are  double  those  for 
the  first  harmonic,  the  actual  dispersion  law  for  the  second 
harmonics  is 

TT2f22    =    gk2/2    +    Tk23/8p. 

Then  the  corrected  second-harmonic  frequency  is 

"l  +  Tk22/pg 


(f2) 


corr 


2  +  Tk  72Pg 


1/2 
f2- 


All  second-harmonic  damping  coefficients  are  plotted  at  their 

corrected  frequencies. 
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The  paddle-wave  results  can  be  fitted  to  the  power-law 
relation 


et  =  cf 

-3      2   -1 
where  c   is  3.2  x  10    sec  cm    (first  harmonics) ,  and  c   is 

-3     2-1 
1.8  x  10    sec  cm    (second  harmonics).   Based  on  the  reasoning 

given  above,  we  view  the  reasonable  agreement  with  Phillips' 
8/3-power  prediction  as  accidental.   Also,  Phillips'  result  rests 
on  a  model  which  is  quite  unlike  the  experiment:  surface  tension 
is  neglected;  the  incident  wave  (already  assumed  small:  ka<<l) 
is  taken  to  vary  little  over  the  time  of  the  decay  calculation; 
the  turbulence  spectrum  is  fixed  on  dimensional  reasoning  based 
on  the  existence  of  an  inertial  subrange.   All  of  these  assump- 
tions are  violated  here. 

The  discrepancy  between  values  of  c  for  the  first  two  har- 
monics suggests  either  an  energy  transfer  between  the  harmonics 
or,  since  the  first  harmonics  are  stronger  than  the  second,  that 
there  is  an  amplitude  dependence  hidden  in  the  damping  coeffi- 
cient B  .   McGoldrick  (1965)  has  shown  that  long-crested  capillary' 
gravity  waves  do  experience  a  resonant  interaction  between  the 
first  and  second  harmonics  near  the  fundamental  wave  number 
K  =  (g/2T)      (2.6  cm    for  pure  water  at  20  C) .   For  measured 
values  of  T,  however,  K  was  far  above  the  experimental  wave  num- 
bers, so  that  resonance  was  highly  unlikely. 


B.   Quadratic  Decay  Hypothesis 
The  alternative  explanation  for  the  discrepancy  between  the 
first-  and  second-harmonic  damping  coefficients  is  explored  by 
assuming  (somewhat  arbitrarily)  a  quadratic  decay  law,  such  as 
might  be  expected  with  turbulent  dissipation,  adding  to  the  linear 
laminar  decay: 

da/dx  =  ~(B£a  +  Yfca2)  (4) 

or 

a (x  )  exp  { -3 „  (x  -  x  )  } 
/     \  o  x.       o  ,  c  . 

a(x)   =  ; : .  5 

Y.a(x  ) 
1  +  -^ — [1  -  exp  {-  B0(x  -  x  )  }] 

The  observed  turbulent  decay  and  eqn.  5  give  y  .   The  results 
are  shown  in  Fig.  6.   Although  there  is  still  some  scatter,  a 
more  consistent  picture  has  emerged.   Aside  from  two  questionable 
experimental  points  (common  to  both  hypotheses) ,  all  points  can 
be  fit  to  the  power-law  relation 

Yt  =  Cf5  (6) 

-5     5-2 
where  C  is  1.7  x  10    sec  cm 

On  the  basis  of  this  evidence,  we  suggest  that  steep  capil- 
lary-gravity waves  in  the  presence  of  underwater  turbulence  obey 
a  dissipative  decay  law  which  is  close  to  quadratic  (i.e., 
eqn.  4)  when 

(i)    the  length  scale  of  the  turbulence  is  much  less  than 
the  wave  length 
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(ii)   the  energy  per  volume  of  the  turbulence  is  the  same  order 
as  the  maximum  energy  per  volume  of  the  waves. 
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List  of  Figures 

Fig.  1.   Schematic  diagram  of  wave  tank  (not  to  scale) .   The  tank 
width  was  25.  cm.   For  wind-driven  waves  a  centrifugal  fan  was 
placed  at  the  paddle  position  and  a  flume  created  by  covering  the 
wave  tank  to  some  distance  past  probe  2. 

Fig.  2.   Upstream  wave  spectrum  obtained  with  paddle  frequency 
3.5  Hz. 

Fig.  3.  Line  spectrum  interpretation  of  wave  spectra  upstream 
(solid  lines)  and  downstream  of  turbulence  (dashed  lines,  dis- 
placed for  clarity).   Paddle  frequency  was  3.1  Hz. 

Fig.  4.   Velocity  component  spectra  2.  cm  below  water  surface  (band 
width  1.0  Hz).   Curve  "a"  with  data  pointsQwas  measured  with  grid- 
generated  turbulence  and  no  waves.   Curve  "b"  with  data  points®, 
was  measured  under  surface  waves  of  paddle  frequency  3.1  Hz,  no 
turbulence.   The  x  indicates  circuit  noise  level. 

Fig.  5.   Turbulent  damping  coefficients,  3  ,  calculated  for  the 
linear  decay  hypothesis,  eqn.  (3),  plotted  as  a  function  of  fre- 
quency.  The  paddle-generated  fundamentals  are  indicated  by O , 
the  second  harmonics  by  x  and  the  wind  generated  wave  measurements 
by  A .   A  "best"  slope  of  2.0  is  drawn  for  comparison. 

Fig.  6.   Turbulent  damping  coefficients,  y  ,  calculated  for  the 
quadratic  decay  hypothesis,  eqn.  (4),  plotted  as  a  function  of 
frequency.   The  paddle-generated  first-harmonics  are  indicated 
by O  ,  the  second  harmonics  by  x,  and  the  wind  generated  waves  by 
A  .   The  straight  line  is  drawn  at  slope  5:1. 
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